Impulsive electron acceleration by Gravitational Waves 
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ABSTRACT 

We investigate the non-linear interaction of a strong Gravitational Wave with the 
plasma during the collapse of a massive magnetized star to form a black hole, or during 



the merging of neutron star binaries (central engine). We found that under certain 
conditions this coupling may result in an efficient energy space diffusion of particles. 
. We suggest that the atmosphere created around the central engine is filled with 3-D 

magnetic neutral sheets (magnetic nulls). We demonstrate that the passage of strong 
pulses of Gravitational Waves through the magnetic neutral sheets accelerates electrons 
ly-j \ to very high energies. Superposition of many such short lived accelerators, embedded 

inside a turbulent plasma, may be the source for the observed impulsive short lived 
bursts. We conclude that in several astrophysical events, gravitational pulses may 



accelerate the tail of the ambient plasma to very high energies and become the driver 
for many types of astrophysical bursts. 



Subject headings: gravitational waves-compact objects: — particle acceleration. 



Introduction 



The interaction of Gravitational Waves (GW) with the plasma and/or the electromagnetic 
waves propagating inside the plasma, has been studied extensively (DeWitt and Breheme 1960 
Cooperstock 1968; Zeldovich 1974; Gerlach 1974; Grishchuk and Polnarev 1976; Denisov 1978 
Macdonald and Thorne 1982; Demianski 1985; Daniel and Tajima 1997; Brodin and Marklund 1999 
Marklund et al. 2000; Brodin et al. 2000; Servin et al. 2000, 2001; Brodin et al. 2001; Moortgat 
and Kuijpers 2003) . All well known approaches for the study of the wave-plasma interaction have 
been used, namely the Vlasov-Maxwell equations (Macedo and Nelson 1982), the MHD equations 
(Papadopoulos and Esposito 1981; Papadopoulos et al. 2001; Moortgat and Kuijpers 2003) and 
the non-linear evolution of charged particles interacting with a monochromatic GW (Varvoglis and 
Papadopoulos 1992). The Vlasov-Maxwell equations and the MHD equations were mainly used to 
investigate the linear coupling of the GW with the normal modes of the ambient plasma, but the 
normal mode analysis is a valid approximation only when the GW is relatively weak and the orbits 
of the charged particles are assumed to remain close to the undisturbed ones. Several studies have 
also explored, using the weak turbulence theory, the non-linear wave-wave interaction of plasma 
waves with the GW (see Brodin et al. (2000)). 
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The strong nonlinear coupling of isolated charged particles with a coherent GW was studied 
using the Hamiltonian formalism (Varvoglis and Papadopoulos 1992; Kleidis et al. 1993, 1995). 
The main conclusion of these studies was that the coupling between GW and an isolated charged 
particle gyrating inside a constant magnetic field can be very strong only if the GW is very intense. 
This type of analysis can treat the full non-linear coupling of the charged particle with the GW 
but looses all the collective phenomena associated with the excitation of waves inside the plasma 
and the back reaction of the plasma onto the GW. 

In this article, we re-investigate the non-linear interaction of an electron with a GW inside a 
magnetic field, using the Hamiltonian formalism. Our study is applicable at the neighborhood of 
the central engine (collapsing massive magnetic star (Fryer et al. 2002; Dimmelmeir et al. 2002; 
Baumgarte and Shapiro 2003)) or during the final stages of the merging of neutron star binaries 
(Ruffert and Janka 1998; Shibata and Uryu 2002). We find that a strong but low frequency (10 
KHz) GW can resonate with ambient electrons only in the neighborhood of magnetic neutral sheets 
and accelerates them to very high energies in milliseconds. Relativistic electrons travel along the 
magnetic field, escaping from the neutral sheet to the super strong magnetic field, and emitting 
synchrotron radiation. We propose that the passage of a GW through numerous localized neutral 
sheets will create spiky sources which collectively produce the highly variable in time. 

In section 2, we analyze the interaction of a GW with a single electron inside a constant 
magnetic field. In section 3 we investigate the diffusion to high energies of a distribution of electrons 
inside the GW. In section 4, we propose a new model for strong coupling of the GW with the 
turbulent plasma at the atmosphere of the central engine and the resulting impulsive synchrotron 
emission and finally, in section 5, we summarize our results. 



2. The Hamiltonian formulation of the GW-particle interaction 

The motion of a charged particle in a curved space and in the presence of a magnetic field is 
described by a Hamiltonian, which, in a system of units m = c = G = 1, is given by 

H(x a ,p a ) = ^g' IU (p^-eA^(p u -eA u ) = ^, a, fi, v = 0, 3. (1) 

g^v _ g^v^ x a^ are foe con travariant components of the metric tensor of the curved space and 
= A^(x a ) are the components of the vector potential of the magnetic field (Misner et al. 
1973). The variables p a are the generalized momenta corresponding to the coordinates x a , and 
their evolution with respect to the proper time r is given by the canonical equations 

dx^_ _ dH_ dpa dH_ 

~drr~^p~ a ' ~dr~~dx^' ^ ' 



We assume a constant magnetic field B = B§e z which is produced by the vector potential 

A = A 1 = A 3 = 0, A 2 = B (x l + c ), c : const., (3) 
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and that a GW propagates in a direction k of angle 6 with respect to the direction of the mag- 
netic field. In that case the nonzero components of the metric tensor are (see Ohanian (1976); 
Papadopoulos and Esposito (1981)) g 00 = 1 and 



11 _ l-asin 2 6* cos ip „22 _ —1 

" — 1+a cost/> " 1+a cost/> (a\ 

33 _ 1-acos 2 6>cos ip 13 _ 31 _ (-a/2) sjnjflcosj/; V > 

y — l+acostf> " y — l+acost/> ' 



where a is the amplitude of the GW and ip = k^x^ = ^(sin^x 1 + cos6>x 3 — x°). The parameter 
v is the relative frequency of the GW, i.e. v = u/fl, where ft = eBo/mc is the Larmor angular 
frequency. We make use of the scaling cBq = 1, thus Vl = 1. 

In the above formalism, the coordinate x 2 is ignorable, so P2 = const.. By setting the constant 
Co in Eq. (3) equal to p2 we get an appropriate gauge that reduces by one degree of freedom the 
Hamiltonian (Eq. (1)), which takes the form 

1 / o 1 — aslcosip 2 1 — acl cos tp 2 2asgcg cos ip x? \ ,_. 

H = - \p - -; —Pi ~ -; —P3 + -PWz ~ t- r , 5 

2 \ 1 — a cos ip 1 — a cos ip 1 — a cos ip 1 + a cos ^> / 

where we use the notation cq = cos and s$ = sin # for brevity. We apply the canonical transfor- 
mation of variables (x , x 1 ,x 3 ,po,pi,p3) — > (x, 9, 0, J) using the generating function 

^(x ,^ 1 ,^ 3 ,/,^, J) = x°I + x x p + z^x 1 + c e x 3 - x°)J. (6) 

The relation between the old and the new variables is given by the equations 

X = x°, I = Po+Pz/ce 

q = x 1 , p = pi - {s e /c e )p 3 (7) 

4> = v(s e x 1 + c e x 3 - x°) J = pz/{c e s e ). 

In the new variables the Hamiltonian (Eq. (5)) takes the form 

1 ( t2 ot,.t o„ ,.t„ l-asicos^ 2 



H = - I I 2 - 2IvJ - 2s e uJp - p 2 - — 1 . (8) 

2 \ 1 — a cos <p 1 + a cos q> J 

Since the variable \ is ignorable, I is a constant of motion and Eq. (8) can be studied as a system 
of two degrees of freedom, where I is a parameter. The variables q and p are associated with the 
gyro-motion. H is of mod(2-7r) with respect to the angle-variable <p and the variable J is related 
linearly with the energy 7 = (1 — v 2 )~ l l 2 of the particles according to the equation 



7 = I - v J. (9) 



The equations of motion are 



T 1— as'i cos l, 

~ s evJ - 1 " ,a P P 



1— acos<j> * " 1+a cos (/> 



4> = -vI-s e vp j = § ( (1+a f os0) . - ( i_l P o 2 s0) . ) sin^ 



(10) 
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where the dot means derivative with respect to the proper time r. Furthermore, Eq.(8) can be 
written as a perturbed Hamiltonian in the usual way, i.e. 

H = H + aH 1 + a 2 H 2 + (11) 

where 

H m = ~(4p 2 + (-l)V) cos m <j>, m > 1 (12) 
are the perturbation terms and 

H = i(/ 2 - 2IuJ - 2s 6 vJ P ) - + Q 2 ) (13) 

is the integrable part of the system that describes the unperturbed helical motion of the particle 
in the flat space. Considering action-angle variables (J±, J 2 , <f>i, 4>2), Eq.(13) takes the form 

H (Ji, <h) = \~ IvJi + ^fjf - J 2 , (14) 

where <f>\ = <f>,J\ = J and 

•h = — f pdq = -(I 2 — 2IvJ + s 2 e v 2 J 2 — 1) , <p2 = arcsin ' 



2vr/ r 2 V ' \V2J2 

Therefore, the unperturbed system is isoenergeticaly non-degenerate for 9 7^ (Arnol'd, et al. 1987) 
and the gyro motion of the particles is represented by trajectories that twist invariant tori with 
angular frequencies u\ = dHo/dJi and u>2 = dE^jdJi- The periodic or quasi-periodic evolution of 
the trajectories depends on whether the rotation number, defined by 

p = — = vl- s 2 e v 2 Ji = v{c 2 e I + S07), (15) 

is rational or irrational, respectively. 

Most of the invariant tori will persist with the presence of the perturbation introduced by 
the GW, if the amplitude is suficiently small, according to the KAM theorem (Arnol'd, et al. 

1987) . The orbits of the particles remain close to the unperturbed ones but their projection on the 
x 1 — x 2 plane is not exactly circular and periodic. Close to the resonant tori, where p is rational, 
the Poincare-Birkhoff theorem applies; a finite number of pairs of stable and unstable periodic 
trajectories survive, producing locally a pendulum like topology in phase space (Sagdeev, et al. 

1988) . 

Since the system is of two degrees of freedom, we can study its evolution by using the Poincare 
sections Ps = {{4>, 7), q = 0, H = 1/2} choosing specific sets of the parameters a, /, v and 9. In the 
numerical calculations, which will follow, we set / = 1. For the unperturbed system (a = 0) the 
sections show invariant curves 7 =const. For o/O some typical examples are shown in Fig.l. 

For small values of a (Fig. la), the invariant curves are perturbed slightly and only close 
to the most significant resonances their deformation becomes noticeable. Increasing further the 
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perturbation parameter a, the width of the resonances increases and homoclinic chaos becomes 
more obvious close to the hyperbolic fixed points (Fig. lb). The existence of invariant curves, which 
confine the resonant regions, guarantees the bounded variation of the particle's energy (A7 = 
0(y/a)) for the chaotic trajectories. 

When the amplitude a exceeds a critical value a c , overlapping of resonances takes place and 
large chaotic regions are generated (Fig.lc) (see also Chirikov (1979)). Particles with initial energy 
7 greater than a critical value 7 C may follow a chaotic orbit which diffuse to regions of higher 
energy, and this will lead them to very high energies in short time scales. For relatively large values 
of a c << a < 1, the islands of regular motion, which survive from the resonance overlapping, 
are gradually destroyed and chaos extends down to relatively low energy particles (Fig. Id). The 
chaotic part of the phase space will be called "the chaotic sea" . 

The dynamics, presented by the Poincare sections in Fig.l, is typical for the majority of 
parameter values. Generally, the critical values a c and j c determine the conditions for possible 
chaotic diffusion. The dynamics of the charged particles shows some exceptional characteristics 
when the frequency of the GW is comparable to the Larmor frequency of the unperturbed motion, 
particularly when 1 < v < 3. For such parameter values, stochastic behavior will appear when 
7 = 1 and for sufficiently large perturbation values large chaotic regions are generated and diffusion, 
even for particles with very low initial energies, will be possible. An example is shown in Fig.2a. 

The evolution of the particles changes character when the direction of propagation of the GW 
is almost parallel to B. In this case, chaos disappears, and the particles undergo large energy 
oscillations. As it is shown in Fig. 2b, a particle, starting even from rest (7 ~ 1), will be driven 
regularly to high energies (7 > 20) and returns back to its initial energy in an almost periodic way. 
In a realistic, non infinite system, several particles may escape from the interaction with the GW 
before returning back to low energies. At 6 = the system is integrable and the energy of the 
particles shows regular slow oscillations with an amplitude proportional to a. 

3. Chaotic diffusion and particle acceleration 

In the previous section, we showed that chaotic diffusion is possible for a > a c and for the 
particles with 7 > j c . Such conditions are necessary but not sufficient for acceleration, since 
islands of regular motion may be present inside the wide chaotic region (see for example Fig. 2). 
We estimate approximately the critical value a c by studying a large number of Poincare sections for 
different values of 9 (Fig.3a). We observe that for small values of v resonance overlap is obtained 
for relatively large amplitudes of the GW. For higher frequency, v > 5, and 9 > 45 the overlapping 
of resonances takes place even for relatively small values of a. The critical value 7 C along the 
frequency axis v is presented for 9 = 45° and for a = 0.05, 0.1, 0.2 in Fig.3b. For large values of u, 
and when the critical value a c is relatively small (~ 0.005), the critical particle energy j c is large 
and increases linearly with v. Namely, for large v, chaotic diffusion takes place only for high energy 
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particles. On the other hand, for v « 2, resonance overlapping takes place for large perturbations 
but the wide chaotic sea formed extends down to the thermal velocity. 

In Fig. 4a the evolution of 7 along a temporarily trapped chaotic orbit (7 < j c ) and an orbit 
which undergoes fast diffusion is shown, using a = 0.02. In Fig. 4b we plot the orbit of a particle 
which on the average is not gaining energy and the average rate of energy gain of 200 particles. 
The diffusion rate of the particles in the energy space is initially fast but for time t > 5000 it starts 
to slow down. The time t is normalized with the gyro period 2-k/VL 

3.1. GW interacting with an ensemble of electrons 

We study next the evolution of an energy distribution N(j, t = 0) of electrons interacting with 
the GW. 

In Fig. 5a we follow the evolution of 3 x 10 4 particles forming initially a cold energy distribution 
N(j,t = 0) ~ 5(7 — 3), where 5 is the Dirac delta function i.e. all particles have the same initial 
energy 7 = 3. A large spread in their energy is achieved in short time scales, and for t = 1000, a 
non-thermal tail extending up to 7 = 100 is formed. 

We repeat the same analysis, assuming that the initial distribution is the tail (v > Vthe, where 
Vthe is the ambient thermal velocity) of a Maxwellian distribution (Fig. 5b). The distribution of the 
high energy particle form a long non-thermal tail analogously to the results reported in Fig. (5a). 

The mean energy diffusion as a function of time is plotted in Fig. 6a for a particular set of 
parameters, and it has the general form 

< 7 >~ t d . (16) 

From a large number of calculations, we find that the energy spread in time follow a normal 
diffusion (d = 0.5) in energy space but as a increases (see Fig. 6b), the interaction becomes super- 
diffusive (d > 0.5) in energy space. This allows electrons to spread fast in energy space and explains 
the efficient coupling between the GW and the plasma. 

4. A model for bursty acceleration 

Our main findings so far are: 

1. The GW can accelerate electrons from the tail of the ambient velocity distribution (i>o(t = 
0) > Vthe) t° very high energies (7 > 10 — 100) for typical values of 0.005 < a < 0.5 and 
5 < uj/VL < 20. Assuming that the the frequency of the GW is around lOKHz, and VL ~ lj/3 
the magnetic fields strength should be around ~ 10~ 4 Gauss. 



- 7- 



2. The acceleration time depends in general on a but is relatively short t acc ~ 100$7 ~ msecs. 
During that time the GW will travel a distance i acc ~ ct acc ~ 10 s cm. 

3. The mean energy diffusion of the electrons interacting with the GW follows the simple scaling 
< 7 >~ t d . 

We propose a new mechanism for efficient particle acceleration around strong and impulsive 
sources of GW using the estimates presented above for the strong interaction of GW with electrons. 
We assume that in the atmosphere of the central engine a turbulent magnetic field will be formed. 
Inside this complex magnetic topology, a distribution of 3-D magnetic neutral sheets (magnetic null 
surfaces) with characteristic length £ will be developed (see Fig. (7) and the relevant studies for the 
solar atmosphere (Bulanov and Sakai 1997; Albright 1999; Longcop and Klapper 2002; Longcop et 
al 2003)). These 3-D surfaces naturally appear and disappear inside driven MHD turbulent plasmas 
and are randomly distributed inside the atmosphere of the central engine. It is well known that 
the neutral sheets will act as localized dissipation regions and will accelerate particles (Nodes et 
al. 2003). The turbulent atmosphere of the central engine will sporadically emit weak X-ray and 
possibly 7-ray bursts for relatively long times before the collapse of the massive magnetized star to 
form a black hole or the merging of neutron star binaries. We suggest that numerous weak bursts 
are present and remain unobserved. 

In this article, we emphasize the role of the GW passing through magnetic neutral sheets and 
claim that the GW will enhance dramatically the acceleration process inside the neutral sheet, 
causing very intense bursts. According to the arguments presented in sections 2 and 3, the GW can 
resonate with the electrons only when the magnetic field is weak (in the vicinity of the magnetic null 
surface). Inside the neutral sheet the synchrotron emission losses are negligible so the acceleration 
is very efficient. Electrons escaping from the neutral sheet to the super strong magnetic field, 
expected in the atmosphere of the central engine, transfer very efficiently their perpendicular to 
the magnetic field energy to synchrotron radiation creating spiky localized bursts. The relativistic 
electrons retain the parallel to the magnetic field energy and travel along the field lines till they 
reach the dense ambient plasma and deposit their energy via collisions, causing other types of longer 
lived bursts (e.g. X-ray and optical flashes) (see Fig. (8a)). 

We suggest thus that the creation of a network of localized accelerators with characteristic 
length £ acc , which are spread in relatively large volumes inside the turbulent atmosphere of the 
central engine, will be responsible for the busty emission. The amplitude of the GW decays as it 
propagates away from the central engine as a ~ 1/r. When a « a c and/or the magnetic topology 
does not allow formation of magnetic neutral sheets the GW stops to interact with electrons and this 
indicates the end of the burst. So the overall duration of the burst is approximately AT ~ AL/c, 
where AL is the length of the characteristic layer where the interaction of the GW with the particles 
is efficient, (see Fig. 8b). 
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4.1. Energy spectrum 

Let us now discuss another very important observational fact. It is well documented that 
the observed spectrum of the synchrotron radiation has a power law shape. This implies that the 
energy distribution causing the synchrotron radiation has also a power law dependence in energy 

7V( 7 ) ~ 7 ~ £ 

were e ~ 1.8 — 2.0 (e.g. see Schaefer et al. (1992) for GRB). These types of energy distribution 
are usually attributed to a shock wave (or a series of shock waves). The formation of a series of 
localized shock waves along the streaming plasma injected from the central engine is considered 
the "standard" model for GRB (Piran 2003; Meszaros 2002) and may be of interest in many other 
astrophysical sources. 

The spontaneous formation of magnetic neutral sheets in the 3-D turbulent atmosphere of the 
central engine is not easy to describe. A number of studies on the statistical properties of the 
magnetic neutral sheets have appeared for the active region of the Sun (Georgoulis et al. 1998; 
Longcop and Noonan 2000; Isliker et al. 2000, 2001; Weatland 2002; Craig and Wheatland 2003). 
All these studies are based on the observational fact that flares follow a specific distribution in 
energy (Crosby et al. 1993). 

An important concept from the study of 3-D MHD turbulence is that a hierarchy of magnetic 
neutral sheets formed, with a distribution of characteristic lengths given by the function 

N(£) ~ r b . (17) 

This scaling follows the energy distribution of flares and the most probable value for b is 3/2. In 
other words, inside a turbulent magnetized plasma, a large number of magnetic neutral surfaces 
with small characteristic lengths, and very few large scale magnetic neutral sheets will be present. 
The exponent, according to the studies reported earlier, for the solar atmosphere is roughly equal 
to b ~ 1 — 2 and is related to the statistics of X-ray bursts (Crosby et al. 1993). The GW will 
accelerate electrons only in the brief period it passes through the magnetic neutral sheet, so the 
acceleration time will also follow a similar scaling law since 

tacc ~ f-acc/c (18) 

We have shown in Fig. 6 that the the mean energy of the accelerated particles increases with 
time according to a simple power law as well. Combining the Eq. (17)-(18) we derive first the 
distribution of t which is a function of I. We have N(t)dt = N(l)dl, or N(t) = N(l)§, which 
yields, when neglecting the constants, 

N(t) ~ t~ b . (19) 

We determine next the distribution of 7 (or, more precisely, of < 7 >), which is a function of 
t. We start again from the relation N(j)d>y = N(t)dt, or N(j) = N{t)§-. We insert Eq. (19), 
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iV(7) = t b ~jt, and we replace t on the right hand side using Eq. (16), 

or, on doing the derivative and rearranging, 

jV( 7 ) = 7 (-6+i-d)/d. (20) 

Using a typical values for d ~ 0.5 (see Fig. 6) and 6 ~ 3/2, we estimate the exponent of the energy 
distribution to be around 2. 



4.2. Energetics 

The energy transferred from the GW to the high energy particles by passing through one 
neutral sheet is 



Ens ~ O n ° ( ^ cc x A£) x im<? (2i) 

where no is the electrons density, rit are the electron at the tail of the local Maxwellian, Al is the 
thickness of the magnetic neutral sheet. The filling factor of neutral sheets inside the atmosphere 
of the central engine is assumed to be /. So the total number of magnetic neutral sheets expected 
to interact with the GW is Al an<1 * ne total energy transferred to the plasma 



tot 



' (AL) 3 



£ 2 x Al 



-n (4c x M) x jmc 2 ) . (22) 
no 



Using typical numbers no = 10 12 cm 3 ,n t /riQ ~ 10 1 ,7 = 100, £ acc ~ 10 8 cm,A£ ~ 10 7 cm / ~ 
10" 1 and assuming that the burst duration is Is which implies that AL ~ 3 x 10 cm, the total 
energy estimated is approximately 10 47 — 10 48 ergs. The total burst duration is around 100s but 
the burst is composed with many short bursts lasting less than a second (£/c). 

We can now list several characteristics of the bursty emission driven by the model proposed 
above: 



• A fraction of the energy carried by the orbital energy of the neutron stars at merger will go 
to the the GW and a portion of this energy will be transferred to the high energy electrons. 

• The topology of the magnetic field varies from event to event, so every burst has its own 
characteristics. 

• The superposition of many small scale localized sources produces a fine time structure on the 
burst. 
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• The superposition of null surfaces with a power law distribution of the acceleration lengths 
will result in a power law energy distribution for the accelerated electrons and an associated 
synchrotron radiation emitted by the relativistic electrons. 

• The decay of the amplitude of the GW and / or the lack of magnetic neutral sheets away from 
the central engine will mark the end of the burst, but not necessarily the end of other types 
of bursts since the cooling of the ambient turbulent plasma has a much longer time scale. 

It is worth drawing the analogy between the elements used to build our model for this busrty 
emission and two significant developments in solar flares and the Earths magnetic tail. (1) Galsgaard 
and Nordlund (1997) studied the dynamic formation of current sheets (current sheets are not always 
associated with neutral sheets) in the solar atmosphere. They found that the irregular motion of 
the magnetic foot-points, caused by the turbulent convection zone, forces the magnetic topologies 
to form a hierarchy of 3-D current sheets inside the complex magnetic topologies above active 
regions. (2) Ambrosiano et al. (1988) studied the passage of low frequency waves (Alfven waves) 
through a 2-D neutral sheet. They discovered that the acceleration of charged particles is greatly 
enhanced from the presence of the Alfven waves. 

We believe that solar flares, the earths magnetotail, and the GW driven bursts share two 
important characteristics, namely the formation of 3-D magnetic null surfaces coupled with the 
passage of waves through these surfaces (Alfven waves for the sun and the earths magnetotail and 
GW for the bursts analyzed in this article). Although the external drivers and the energetics of 
the bursts produced are clearly very different. 



5. Summary 

In this article we have study an efficient mechanism to transfer the energy released from 
a central engine to the turbulent plasma in its atmosphere. We have analyzed the non-linear 
interaction of GW with electrons. We showed that electrons in the tail of the ambient plasma are 
accelerated efficiently, reaching energies up to hundreds of MeV in less than a second, when (1) the 
amplitude of the GW is above a typical value 0.005 — 0.5 and (2) the magnetic field is relatively 
weak Bq ~ 10 -4 Gauss, for a GW with typical frequency l-10KHz. 

On the basis of these findings, we propose that pulsed GW emitted from the central engine 
will interact with the ambient plasma in the vicinity of the magnetic neutral sheets formed natu- 
rally inside externally driven turbulent MHD plasmas. Magnetic neutral sheets have characteristic 
lengths I ~ 10 7 - 10 s cm and are short lived 3-D surfaces. Although these structures are efficient 
accelerators, we are emphasizing in this article only the role of the GW passing through these sur- 
faces since we focus our attention on the very strong and bursty sources. The GW passing through 
the neutral sheets will accelerate electrons to very high energies. Relativistic electrons escape from 
the magnetic neutral sheets radiating synchrotron emission as soon as they reach the very strong 
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magnetic fields. 

The collective emission of thousands of short lived (less than a second) synchrotron pulses, 
created during the passage of the GW through a relatively large volume. The relativistic electrons 
loosing most of their perpendicular to the magnetic field energy to synchrotron radiation retain 
the parallel energy and heat the ambient plasma emitting other types of longer lived bursts, e.g. 
X-rays and optical flushes. 

We have also shown that if the acceleration length follows a simple scaling law (N(£ acc ) ~ (-ace) 
and since, as shown, < 7 >~ t^ cc , the energy distribution also follows a power law scaling (N(< 
7 >) ~< 7 >( _b + 1_c 0/ (i ) ) which for reasonable values of b ~ 3/2 and d ~ 1/2 agrees remarkably 
well with the energy distributions inferred from the observations. 

We are suggesting that numerous burst are still present before and after the passage of the 
GW since the magnetic neutral sheet act as efficient particle accelerators. The GW passage simply 
enhances the acceleration process and makes many short lived spikes to become visible. We pro- 
pose that a more thorough analysis of the statistical properties of the bursty emission, with main 
emphasis on the low energy part of the frequency distribution, will give important insight on the 
processes mentioned here. 

A detailed model for the interaction of GW with turbulent MHD plasma is currently under 
study, and we hope to develop an even more efficient energy transfer from the GW to the plasma e.g 
by triggering the interaction (percolation) of many null sheets during the passage of the GW. We 
hope that this may lead us to an alternative scenario for the still unresolved questions related with 
the acceleration mechanism in the atmosphere of the central engines and the physical processes 
behind the X-ray and GRB. 
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Fig. 4- 
2.2,0(0) 



a) The evolution of 7 along a trapped in a magnetic island chaotic orbit for 7(0) 



7T and along a diffusive one for 7(0) = 2.6, 0(0) = (a = 0.02, (9 = 45 c 



5) b) The 



evolution of 7 along a strongly chaotic orbit (dotted line) and its average value (solid line) along 
200 trajectories starting with 7(0) = 2.0 and a randomly selected 0(0) (a = 0.1,8 = 45°, v = 5). 
The time is normalized with the gyro-period (2tt/S1). 



-16- 



10000 



1000 




100 - 



10000 



1000 



100 - 



10 - 




Fig. 5. — The evolution of an energy distribution, a) The initial distribution (dotted line) consists 
of 3 x 10 particles having y(t = 0) = 3. b) The initial distributions is Maxwellian, as it is shown 
by the dotted curve. Only particles in the tail of the Maxwellian with 7 > 7 C will be accelerated. 
The parameters used in both studies are a = 0.5, v = 20 and = 30°. 
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Fig. 6. — (a) We plot the mean energy diffusion in time. We observe that log(< 7 >) is linearly 
related to log(t) for t > t*. The slope d for a = 0.7 is d ~ 0.6. (b) The slope d is plotted as a 
function of the GW amplitude. 



-17- 




Fig. 7. — A schematic representation of the thin three dimensional magnetic null sheets appearing 
spontaneously and fill densely a driven turbulent magnetized plasma 



18 - 




(a) (b) 

Fig. 8. — (a) The propagation of a GW through a magnetic neutral sheet accelerates electrons very 
efficiently. Relativistic electrons stream away from the accelerator and emit a pulse of synchrotron 
radiation when they reach the super strong magnetic fields. The dashed line represents the magnetic 
field null surface and l acc is the acceleration length, (b) A collection of magnetic neutral sheets is 
formed inside the turbulent atmosphere (region C) of the central engine (region A). A GW pulse 
propagating away from the central engine and passing through the region C will form numerous 
7-ray spikes by accelerating particles near the magnetic null surfaces. The superposition of these 
spikes form a short lived burst. The total burst duration is approximately 100s (AT ~ AL/c) but 
it is composed by many short spikes lasting less than a second (i/c). The GW pulse will become 
very weak and the density of the magnetic null surfaces will drop dramatically in the region D, and 
this will mark the end of the burst. 



